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ABSTRACT 

SCF-Xo-SW  molecular-orbital  calculations  have  been  carried  out  for  several 
configurations  of  hydrogenated  silicon  clusters  in  order  to  determine  the  con¬ 
tribution  of  three-center  bonding  to  the  electronic  structure  of  hydrogenated 
amorphous  silicon.  Three-center  bonding  of  dissociated  molecular  hydrogen  is 
shown  to  stabilize  Si-Si  bonds  over  a  wide  range  of  Si-Si  distances.  It  can 
be  concluded  that  hydrogenation  can  in  principle,  saturate  all  strained  as  well 
as  dangling  bonds  in  a-Si.  The  results  further  indicate  that  a  single  hydrogen 
three-center  bond  is  unlikely  in  a-Si:H  alloys. 
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I.  Introduction 

The  discovery  that  a-$i:H  alloys  can  be  doped  either  n-type  or  p-type^  and 

2 

can  be  used  as  the  active  material  in  high-efficiency  solar  cells  has  stimulated 

a  great  deal  of  interest  in  the  electronic  structure  of  these  semiconductors. 

3  4 

Despite  this  enormous  interest,  however,  many  unanswered  questions  remain.  * 

It  is  clear  that  pure  a-Si  contains  large  defect  concentrations  including 
both  strained  and  dangling  bonds,  and  that  both  of  these  are  compensated  by 

5 

the  addition  of  hydrogen.  Sol  et  aK  have  shown  that  a-Si  deposited  via  the 
thermal  decomposition  of  SiH^  contains  the  order  of  10  unpaired  spins  /cm  , 
but  that  these  are  eliminated  by  introducing  approximately  the  same  concentration 
of  hydrogen.  Furthermore,  they  were  able  to  incorporate  over  100  times  more 
hydrogen  into  such  films  after  removal  of  the  unpaired  spins.  This  is  strong 
evidence  that  spinless  defects  are  present  in  a-Si  films.  Further  evidence  comes 
from  effusion  experiments  which  indicate  that  about  100  times  as  much  hydrogen 
is  given  off  than  unpaired  spins  are  created.** 

There  are  several  possible  explanations  of  the  spinless  defects.  The  most 
straightforward  of  these  is  simply  an  array  of  strained  bonds,  either  with 

stretched  distances  or  distorted  bond  angles.  Other  possibilities  include 

7  8  9 

charged  pairs  of  dangling  bonds,  *  doubly  coordinated  silicon  atoms,  or 

three-center  bonds  with  a  negative  correlation  energy. These  proposed 
defects  have  very  different  electronic  structures,  which  might  allow  for  an 
experimental  determination  of  their  concentrations.  Strained  bonds  would 
be  expected  to  yield  valence  and  conduction  band  tails,  doubly  coordinated 
silicon  centers  should  provide  two  different  EPR  lines  upon  n  and  p  doping, 
and  negatively  correlated  defects  should  tend  to  pin  the  Fermi  energy.  How¬ 
ever,  at  present,  no  experiments  have  been  carried  out  which  can  unambiguously 
eliminate  any  of  the  proposed  spinless  defects. 


In  a  previous  paper,  the  results  of  self-consistent-field  Xa  scattered- 
wave  (SCF-Xo-SW)  molecular-orbital  calculations  on  silane  molecules  and 
silicon-hydrogen  clusters  demonstrated  that  such  clusters  serve  as  models  for 
the  local  electronic  structure  of  a-Si:H  alloys.  However,  the  clusters  in® 
vestl gated  only  contained  terminal  hydrogen  atoms.  In  this  paper,  we  present 
the  results  of  SCF-Xa-SW  calculations  on  clusters  containing  three-center  bonds. 

The  results  of  these  calculations  are  presented  in  Section  II.  In  Section  III, 
we  discuss  the  insights  that  these  results  yield  with  regard  to  the  electronic 
structure  of  a-Si:H  alloys. 

II.  Results 

Before  an  evaluation  of  three-center  bonding  in  a-Si:H  is  possible,  we 
must  model  the  Si-Si  bond  in  these  alloys.  In  order  to  do  this,  we  use  the 
cluster  SigXg  where  X  represents  a  saturator;  Si2Hg  is  an  disilane-like 
molecule  saturated  by  six  hydrogen  atoms  at  a  Si -H  distance  equal  to  the 

O 

crystalline  Si -Si  mea rest-neighbor  separation  of  2.35A.  If  such  a  cluster  is 
to  serve  as  a  good  model  for  the  Si-Si  bond  in  a-Si:H,  the  molecular-orbital 
energy  levels  should  be  insensitive  to  the  conformation  of  the  saturators  yet 
highly  sensitive  to  the  Si-Si  distance.  In  order  to  test  the  model,  we  calculated 
the  molecular-orbital  energy  levels  for  both  the  eclipsed  and  staggered  con¬ 
formations  of  the  saturators  with  the  Si -Si  distance  being  held  at  2.35A,  and 
then  recalculated  the  energy  levels  using  a  staggered  cluster  with  a  varying 
Si-Si  distance.  The  results  are  shown  in  Fig.  1.  It  is  clear  from  a  comparison 
of  Fig.  1(a)  and  1(b)  that  the  energy  levels  are  completely  Insensitive  to  the 
saturator  conformation,  while  Fig.  1(c)  strikingly  shows  their  extreme  sensitivity  to 
a  changing  S1-S1  bond  length.  A  separate  calculation  demonstrates  that  an  in¬ 
crease  of  the  Si-saturator  bond  length  parallels  the  behavior  of  an  increasing 

Si -Si  bond  length.  Thus,  these  clusters  serve  as  good  models  for  the  Si -Si  bond 


in  tetrahedral ly  coordinated  a-Si  alloys. 

Since  we  are  modeling  only  the  Si-Si  bond  rather  than  the  local  atomic 
configuration,  we  can  restrict  our  attention  to  those  molecular  orbitals  that 
have  appreciable  electronic  density  in  the  region  between  the  silicon  atoms. 

These  are  indicated  in  Fig.  1  by  solid  lines  (The  molecular  orbitals  corresponding 
to  the  dotted  lines  are  primarily  located  near  the  saturators.}.  Four  molecular 
orbitals  are  thus  of  interest.  The  lowest-lying  is  gerade  and  is  composed  pre¬ 
dominantly  of  silicon  3s  orbitals.  Since  the  energy  of  this  orbital  is  in¬ 
dependent  of  the  Si-Si  bond  length,  it  is  nonbonding.  The  next- lowest- lying 
molecular  orbital  is  ungerade,  but  is  also  composed  predominantly  of  silicon 
3s  orbitals  and  is  nonbonding.  The  other  two  are  the  Si-Si  bonding  and  anti¬ 
bonding  orbitals.  It  is  clear  from  Fig.  1  that  for  Si-Si  separations  larger 

O 

than  the  crystalline  value  of  2.35A  these  can  have  energies  within  the  band 

O 

gap.  When  the  Si-Si  distance  is  less  than  3.1A,  the  bonding  orbital  is  occupied 
and  the  antibonding  orbital  is  empty,  resulting  in  weak  bonds.  Alternatively, 

o 

for  Si -Si  separations  greater  than  3.1 A,  the  anti  bonding  orbital  is  the  occupied 
one,  yielding  dangling  bonds. 

It  is  generally  assumed  that  hydrogenation  of  pure  a-Si  results  in  the 
saturation  of  the  dangling  Si-Si  bonds  via  the  dissociation  of  molecular  hydrogen 
to  form  two  Si-H  bonds.  This  process  is  sketched  in  Fig.  2.  The  H-H  bonding 
orbital  mixes  with  the  Si-Si  gerade  orbitals  to  form  typel  and  type 3 
molecular  orbitals,  while  the  H-H  antibonding  orbital  mixes  with  the  Si-Si  un¬ 
gerade  orbitals  to  form  type  2  and  type  4  molecular  orbitals.  The  relative 
energies  of  these  molecular  orbitals  depends  on  the  fraction  of  nonbonding 
character  In  the  original  orbitals,  In  turn  a  sensitive  function  of  the  inter¬ 
atomic  distances.  Fig.  3  shows  the  energies  of  the  molecular  orbitals  which 
arise  from  the  dissociation  of  a  hydrogen  molecule  to  saturate  the  dangling  bonds 
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of  a  pair  of  silicon  atoms  3.2A  apart.  The  H-H  bonding  orbital  mixes  with 
the  unoccupied  Si-Si  bonding  orbital,  yielding  a  low-lying  type  1  molecular 
orbital  which  tends  to  widen  the  valence  band  of  the  a-Si.  The  second- lowest- 
lying  molecular  orbital  is  type  4  and  results  from  the  mixing  of  the  ungerade 
Si -Si  nonbonding  and  the  H-H  anti bonding  orbitals.  The  gerade  Si -Si  nonbonding 
orbital  and  the  H-H  bonding  orbital  mix  to  produce  a  type  3  molecular  orbital. 

Finally,  the  Si-Si  antibonding  orbital  mixes  with  the  H-H  antibonding  orbital, 
yielding  a  type 2  molecular  orbital.  The  overall  result  is  the  occupation 
of  two  H-H  bonding  and  two  H-H  anti  bonding  orbitals,  yielding  a  nonbonding  H-H 
interaction,  but  simultaneously  inducing  a  strongly  bonding  Si-H  interar+ion. 

Contour  plots  of  the  four  molecular  orbitals  are  shown  in  Figs.  4  through  7.  It  is  clea 
that  these  results  are  consistent  with  the  idea  of  hydrogen  saturating  the 
dangling  bonds  of  pure  a-Si. 

O 

For  Si-Si  separations  less  than  3.1A,  it  might  appear  that  there  is 
insufficient  room  to  insert  a  hydrogen  molecule  between  the  Si-Si  bond.  However, 
since  the  hydrogen  orbitals  are  of  s  character,  they  are  insensitive  to  the 
angular  position  of  the  atoms.  In  fact,  by  rotating  the  hydrogen  atoms  off 
the  Si -Si  axis,  thus  varying  the  Si -H-H  bond  angle,  the  equilibrium  Si-H  and 
H-H  bond  lengths  can  be  preserved  over  a  wide  range  of  Si -Si  distances.  The 

O 

molecular-orbital  levels  for  fixed  Si-H  separations  of  1.35A  and  H-H  separations 

o  o  o 

of  0.8A  Is  shown  In  Fig.  8  as  the  Si-Si  distance  varies  from  3.4A  to  2.5A.  At 
all  SI-SI  separations,  the  introduction  of  hydrogen  stabilizes  the  orbitals  that 
would  lie  In  the  band  gap  in  pure  a-Si.  Thus,  hydrogen  can  compensate  stretched 
as  well  as  dangling  bonds  In  a-SI. 

A  contour  plot  of  a  stabilized  bonding  orbital  for  a  Si -Si  separation 

O 

of  3.2A  Is  given  In  Fig.  9.  The  stabilization  results  from  Increased  bonding 
between  the  silicon  atoms  and  the  more  distant  hydrogen  atoms,  i.e.  Increased 
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three-center-bonding  activity.  Note  that  when  the  hydrogen  atoms  are  not 
located  along  the  Si-Si  axis,  at  least  two  completely  equivalent  positions 
exist  for  them,  and  we  should  expect  rapid  oscillations  between  these  equivalent 
conformations. 

Fisch  and  Licciardello^  suggested  that  a  single  hydrogen  atom  could 
bridge  between  two  silicon  atoms  in  a-Si:H,  and  they  identified  two  occupied 
molecular  orbitals  with  appreciable  electron  density  in  the  region  between  the 
two  silicons.  In  Fig.  10,  we  show  the  molecular  orbital  levels  which  would 

o 

result  from  a  single  hydrogen  atom  between  two  silicon  atoms  3.2A  apart.  Clearly, 
four  rather  than  two  molecular  orbitals  are  important.  The  lowest-lying  is 
bonding  between  both  silicons  and  the  hydrogen  atom  via  the  Si  3s  orbitals. 

The  second-lowest-lying  is  bonding  between  the  hydrogen  atom  and  the  nearest 
silicon  atom  but  weakly  antibonding  between  the  two  silicon  atoms  (When  the 
hydrogen  atom  is  half-way  between  the  two  silicons,  this  molecular  orbital  is 
Si -H  nonbonding).  The  third-lowest-lying  molecular  orbital  is  bonding  between 
both  silicon  atoms  and  the  hydrogen  atom  via  the  Si  3p  orbitals.  Finally,  the 
highest-lying  occupied  molecular  orbital  is  antibonding  between  the  silicon 
atoms  and  nonbonding  between  the  hydrogen  atom  and  both  silicon  atoms  when  the 
hydrogen  is  half-way  between  the  two  silicons.  However,  if  the  hydrogen  atom 
is  displaced  toward  either  silicon  atom,  this  molecular  orbital  remains  anti¬ 
bonding  between  the  silicons  but  becomes  bonding  between  the  hydrogen  atom  and 
the  nearer  silicon  atom.  This  additional  bonding  interaction  stabilizes  the 
asymmetric  displacement  of  the  hydrogen. 

III.  Discussion 

The  molecular-orbital  stabilization  described  at  the  end  of  Section  II 
has  been  interpreted^  as  leading  to  a  negative  effective  correlation  energy 
for  the  S1-H-S1  defect  center  in  a  Si : H  However,  the  existence  of  a  negative 
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effective  correlation  energy  requires  that  the  energy  to  place  a  second 
electron  in  the  molecular  orbital  be  lower  than  that  of  the  first.  Our 
results  suggest  that  even  if  the  molecular  orbital  is  singly  occupied  the 
hydrogen  atom  is  displaced  from  the  center  because  of  the  increased  bonding 
interactions.  This  stabilization  has  nothing  to  do  with  a  negative  effective 
correlation  energy.  Although  doubly  occupying  the  molecular  orbital  may  result 
in  an  additional  stabilization  from  an  increase  in  the  "hydrogen  bonding"  con¬ 
tribution,  it  would  appear  unlikely  that  this  could  overcome  the  extra  repulsion 
from  the  presence  of  two  electrons  in  the  same  orbital.  In  any  event,  the 
much  greater  stabilization  from  the  presence  of  two  hydrogen  atoms  between 
the  separated  silicon  atoms  should  lead  to  a  negligibly  small  density  of  the 
single-hydrogen  centers,  particularly  in  view  of  the  fact  that  film  growth  in 
a-Si:H  almost  certainly  takes  place  primarily  via  SiH2  reactions.  In  fact, 
the  concerted  cleavage  of  hydrogen  molecules  to  form  metal -hydrogen  bonds,  as 

exemplified  by  the  process  sketched  in  Fig.  3,  is  common  in  many  homogeneous 

13 

catalytic  processes.  All  of  these  homogeneous  processes  require  that  the 

metal  expand  its  valence  shell.  It  is  well  known  that  silicon  reacts  almost 

exclusively  by  expanding  its  valence  shell  and  passing  through  a 

14 

transition  state  of  the  type  shown  in  Fig.  3 .  This,  coupled  with  the  fact 
that  the  concerted  addition  of  hydrogen  molecules  can  stabilize  Si-Si  bonds  over 
a  wide  range  of  Si -SI  separations,  supports  the  conjecture  that  three-center 
bonds  In  a-S1:H  alloys  almost  exclusively  involve  hydrogen  molecules  rather  than 
individual  atoms. 

The  possibility  of  significant  densities  of  three-center  bonds  of  the 
type  described  here  has  not  heretofore  been  considered  in  a-S1:H  alloys.15  A 
major  consequence  out  esults  is  that  stretched  bonds  can  be  compensated  by 
hydrogen  removing  states  from  the  valence  and  conduction  band  tails.  This 
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is  consistent  with  the  sharp  absorption  edges  observed  in  a-Si :H  alloys.16  It 

further  suggest  that  mechanisms  other  than  the  breaking  of  stretched  bonds1 ^ 

1 8 

be  sought  to  account  for  the  observed  photo-induced  changes  which  characterize 
these  materials. 

Three-center  bonds  could  well  be  the  origin  of  one  of  the  peaks  observed 

19 

in  hydrogen  effusion  studies.  The  residual  bonding  component  after  release 

of  a  hydrogen  molecule  from  the  centers  considered  here  suggest  that  effusion 

can  occur  at  a  lower  activation  energy  than  from  either  Si-H^  or  Si -H  centers. 

We  have  not  calculated  the  vibrational  structure  of  the  three-center 

bonds  proposed  herein.  There  have  been  problems  with  unambiguous  identification 

of  the  2000cm’1  infrared  absorption  with  Si -H  centers  and  the  2090cm’1  absorption 
20 

with  Si -Hg  centers.  It  is  important  to  consider  the  additional  infrared  absorption 
from  three-center  bonds  before  such  problems  can  be  fully  resolved. 
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Figure  1:  SCF-Xa-SW  molecular-orbital  energy  levels  of  an  Si 2(sat)g  cluster 
having  the  following  configurations: 

o 

A.  Staggered  conformation  with  Si -Si  distance  =  2.35A 

B.  Eclipsed  conformation  with  Si -Si  distance  =  2.35A 

o 

C.  Staggered  conformation  with  Si -Si  distance  =  3.2A 

The  highest  occupied  molecular  orbitals  are  indicated  by  arrows, 
and  orbitals  with  appreciable  density  between  Si  atoms  are  shown 
as  solid  lines. 

Figure  2:  Bond  types  for  the  formation  of  a  SigfsatL^  complex.  Vertical 
lines  denote  nodes  and  horizontal  lines  indicate  bonding. 

Figure  3:  SCF-Xa-SW  molecular-orbital  energy  levels  of  an  S^Csat)^  complex 

o 

as  a  function  of  H-H  distance,  with  the  Si-Si  distance  held  at  3.2A. 

A.  Staggered  conformation  of  Si 2( sat)g  cluster  shown  for  comparison 
(same  as  Figure  1C). 

B.  H-H  distance  =  0.53A 

C.  H-H  distance  =  0.80A 

D.  H-H  distance  =  1.00A 

The  highest  occupied  molecular  orbitals  are  indicated  by  arrows, 
and  orbitals  with  appreciable  density  between  Si  atoms  are  shown 
as  solid  lines. 

Figure  4:  Contour  map  for  the  Si 2( sat )^H2  molecular-orbital  wave  function  of 
lowest  energy  in  Figure  3C.  This  orbital  is  bonding  between  the 
hydrogen  atoms. 

Figure  5:  Contour  map  for  the  SigfsatJgHg  molecular-orbital  wave  function  of 
second-lowest  energy  in  Figure  3C.  The  solid  and  dashed  contours 
represent  positive  and  negative  vaues,  respectively,  of  the  wave 
function.  This  orbital  is  bonding  between  the  hydrogen  atoms. 


Figure  6:  Contour  map  for  the  Si2(sat)gH2  molecular-orbital  wave  function  of 
third-lowest  energy  in  Figure  3C.  The  solid  and  dashed  contours 
represent  positive  and  negative  values,  respectively,  of  the  wave 
function.  This  orbital  is  anti  bonding  between  the  hydrogen  atoms. 
Figure  7:  Contour  map  for  the  Si2(sat)gH2  molecular-orbital  wave  function  of 
fourth-1 owest  energy  in  Figure  3C.  The  solid  and  dashed  contours 
represent  positive  and  negative  values,  respectively,  of  the  wave 
function.  This  orbital  is  antibonding  between  the  hydrogen  atoms. 
Figure  8:  SCF-Xa-SW  molecular-orbital  energy  levels  of  an  Si'2(sat)gH2  complex 
as  a  function  of  Si -Si  distance  and  H2  molecular  orientation,  with 

o  o 

the  H-H  and  Si-H  distances  held  at  0.8A  and  1.35A,  respectively. 

O 

A.  Si-Si  distance  «  3.4A 

o 

B.  Si-Si  distance  =  3.2A 

o 

C.  Si-Si  distance  *  2.5A 

The  highest  occupied  molecular  orbitals  are  indicated  by  arrows, 
and  orbitals  with  appreciable  density  between  Si  atoms  are  shown 
as  solid  lines. 

Figure  9:  SCF-Xa-SW  molecular-orbital  wave  function  of  lowest  energy  in 

Figure  8B.  The  solid  and  dashed  contours  represent  positive  and 
negative  values,  respectively,  of  the  wave  function.  This  orbital 
is  bonding  between  the  hydrogen  atoms. 

Figure  10:  SCF-Xa-SW  molecular-orbital  energy  levels  of  an  S12(sat)gH  complex 
as  a  function  of  the  position  of  the  H  atom  between  the  Si  atoms, 
with  the  Si -Si  distance  held  at  3.2A. 

A.  SI-H  distance  -  1.6A 

B.  Si-H  distance  -  1.34A 

C.  Si-H  distance  ■  1.2A 

D.  Si-H  distance  -  1.07A 


The  highest  occupied  molecular  orbitals  are  indicated  by  arrows, 
and  orbitals  with  appreciable  density  between  Si  atoms  are  shown 
as  solid  lines. 
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